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Abstract 

This article presents an analysis of about 29,000 measurements of gamma radiation associated with the decay of radon in a sealed 
container at the Geological Survey of Israel (GSI) Laboratory in Jerusalem between 28 January 2007 and 10 May 2010. These 
measurements exhibit strong variations in time of year and time of day, which may be due in part to environmental influences. 
However, time-series analysis reveals a number of periodicities, including two at approximately 11.2 year -1 and 12.5 year . We 
have previously found these oscillations in nuclear-decay data acquired at the Brookhaven National Laboratory (BNL) and at the 
Physikalisch-Technische Bundesanstalt (PTB), and we have suggested that these oscillations are attributable to some form of solar 
radiation that has its origin in the deep solar interior. A curious property of the GSI data is that the annual oscillation is much 
stronger in daytime data than in nighttime data, but the opposite is true for all other oscillations. This may be a systematic effect 
but, if it is not, this property should help narrow the theoretical options for the mechanism responsible for decay-rate variability. 
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1. Introduction 

Radon ( 222 Rn) is a radioactive inert gas formed from the dis- 
integration of 226 Ra as part of the 238 U decay chain. It occurs 
with varying concentrations in the geological environment, and 
is widely used for tracking temporally varying geological pro- 
cesses [1]. The nature of the physical processes driving the 
temporal patterns observed in 222 Rn time series is not clear, 
particularly the extent to which environmental parameters, such 
as temperature, atmospheric pressure, and humidity, can influ- 
ence the 222 Rn level J2j,|3[]. To help understand these processes, 
Steinitz et al. have carried out experiments that track the vari- 
ation of alpha and gamma radiation arising from a source of 
222 Rn in a closed volume of air 10]. In this article, we analyze a 
sequence of 28,733 measurements of gamma radiation acquired 
in this experiment over the time interval 28 January 2007 to 1 
May 2010. The experimental setup is described briefly in Ap- 
pendix A, and environmental drivers of radon production are 
discussed in Appendix B. 

In Section 2, we present the gross features of the data se- 
quence, which exhibits a strong dependence on time of day and 
time of year. We also show similar displays of the laboratory 
temperature and of the solar elevation, finding a stronger corre- 
lation of the gamma-detector measurements with the latter than 
with the former. In Section 3 we present an analysis of the 
measurements in terms of frequency and time of day, finding 
that most of the modulation is associated with nighttime mea- 
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surements. For this reason, we also present separate power- 
spectrum analyses of the daytime and nighttime measurements. 
In Section 4, we compare these results with those of our pre- 
vious time-series analyses of BNL and PTB data JUS 01, and 
we comment on a recent study of the decay rate of I37 Cs in an 
experiment carried out at the Gran Sasso Laboratory [g). Our 
results are discussed in Section 5. 



2. Overview of the Sequence of Gamma-detector measure- 
ments 

The mean daily count rate of photons detected by the gamma 
detector is 2.3610 x 10 5 . These counts are due primarily to the 
beta decay of 214 Pb and of 214 Bi, but there are other sources of 
gamma rays in the 222 Rn decay chain, and one cannot rule out 
that some counts may be due to the decay of the other 222 Rn 
daughters. Figure[T]shows a plot of the daily count, normalized 
to mean value unity, as a function of time. Since the early values 
are anomalous, we exclude them from further consideration. 

We next remove fluctuations on the time scale of one day 
by forming 21 -point running means. Figure 2 shows the per- 
centage departure from the mean value (unity). The dominant 
feature is an annual modulation, with a peak at mid-summer. 
However, it is notable that there is also a minor peak near mid- 
winter. 

Due to the experimental setup and location, it is reasonable to 
suspect that an annual modulation may have an environmental 
origin. We therefore show comparable plots derived from the 
absolute temperature (mean value 291.87 K) in Figure 3; pres- 



Accepted far publication in Astroparticle Physics 



March 3, 2013 



1.15r 
1.1 - 




o 

0.85 
0.8- 

0.75.' ' ' ' ' ' ' 

2007 2007.5 2008 2008.5 2009 2009.5 2010 2010.5 
Date 

Figure 1 : Gamma measurements, normalized to mean value unity, as a function 
of time. 
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Figure 2: Fluctuation (in percent) of the 21 -point running mean of the (abbre- 
viated) gamma measurements. 



sure (mean value 923.78 milli-bar) in Figure 4; and line voltage 
(mean value 13.403 Volts) in Figure 5. 

Pressure shows little annual variation, but temperature and 
voltage both show annual modulations with amplitudes and 
phases comparable with that of the gamma-detector measure- 
ments, except that the voltage curve is in anti-phase. Table [2] 
lists the amplitude of the sine-wave component of each of these 
wave-forms, together with the phase-of-year of the maximum 
and the phase-of-year of the minimum for each sine-wave fit. 
We see that the amplitudes of the temperature and voltage varia- 
tions are comparable with that of the gamma-detector measure- 
ments. However, even allowing for an inversion of the voltage 
curve, the phases differ by over one month. Such a large phase 
shift would be difficult to understand if the annual variation of 
the gamma-detector measurements were caused by the annual 
variations of either the temperature or the voltage. 
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Figure 3: Fluctuation (in percent) of the 21 -point running mean of the absolute 
temperature. 

8i 1 1 1 1 1 1 1 

6- 



4- 




a 

S -2- 



-4- 
-6- 



2007 2007.5 2008 2008.5 2009 2009.5 2010 2010.5 

Date 

Figure 4: Fluctuation (in percent) of the 21-point running mean of the pressure. 



These phase relationships are also evident in Figure 6 and 7. 
Figure 6 combines the waveforms of the gamma-detector mea- 
surements and of the temperature. Figure 7 combines the wave- 
form of the gamma-detector measurements and the inverted 
waveform of the voltage. It is notable that the temperature and 
voltage curves lag the gamma curve. 

The variation of the gamma-detector measurements as a 
function of both date and time of day is shown in Figure 8. 
We see that there is a strong dependence on time of day, as well 
as on time of year. The variation of temperature as a function of 
date and time of day is shown in Figure 9. The corresponding 
plot for solar elevation is shown in Figure 10. The solar eleva- 
tion has its maximum value close to mid-year and mid-day, as 
expected. The peaks in the gamma-detector measurements oc- 
cur close to mid-year, but a little later than mid-day. The peaks 
in the temperature measurements occur slightly later than mid- 
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Figure 5: Fluctuation (in percent) of the 21 -point running mean of the voltage. 



Table 1: Amplitudes and Phases of Gamma, Temperature, Pressure, and Volt- 
age^ 





Amplitude 


Phase of Maximum 


Phase of Minimum 


Gamma 


0.0208 


0.4664 


0.9664 


Temperature 


0.0287 


0.5660 


0.0660 


Pressure 


0.0036 


0.9895 


0.4895 


Voltage 


0.0196 


0.0770 


0.5770 



year and slightly later than mid-day. 

These results suggest that the variation in the gamma- 
detector measurements is related to environmental conditions, 
but not in exactly the same way as is temperature. 

3. Power Spectrum Analyses 

We have carried out a power-spectrum analysis of the 
gamma-detector measurements when broken down by hour of 
day. For the frequency band 0-5 year 1 , the result is shown in 
Figure 1 1 . The strongest feature corresponds to an annual os - 
dilation centered on mid-day. The next strongest feature (fre- 
quency 2 year -1 ) is a semi-annual oscillation. This is mainly a 
night-time feature, but it has a peak at 17 hours. 

Figure 12 covers the frequency band 5-10 year -1 . The oscil- 
lations in this band all occur at nighttime. The strongest feature 
is at frequency 8.9 year -1 . 

Figure 13 covers the frequency band 10-15 year -1 . Here also, 
the oscillations all occur at nighttime. The strongest feature is at 
frequency 12.5 year -1 , but there is also a feature at 1 1.3 year -1 
and a lesser feature at 1 1.9 year -1 . 

It is clear from the preceding results that the power-spectrum 
properties of the gamma-detector measurements will be quite 
different for daytime measurements and nighttime measure- 
ments. We therefore examine separately data acquired over the 
time interval 6 to 18 hours (daytime), and from 18 hours to 6 
hours (nighttime). 

We show in Figure 14 the power spectra derived from day- 
time and nighttime measurements over the frequency range 0-5 
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Figure 6: Fluctuation (in percent) of the 21 -point running mean of the absolute 
temperature (red) and of the gamma measurements (green). 
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Figure 7: Fluctuation (in percent) of the 21 -point running mean of the voltage 
(inverted, red) and of the gamma measurements (green). 

year -1 . We see that there is a very strong peak at 1 year -1 in 
the daytime data, and a strong peak at 2 year -1 in the nighttime 
data, as we would expect from Figure 1 1 . 

We see from Figure 15 that, for the frequency range 5-20 
year -1 , the nighttime data lead to significantly stronger oscil- 
lations than do the daytime data. However, the fine structure 
in these power spectra seems unrealistic, since the duration of 
the time series is only about 3 years. For this reason, we have 
formed a 41 -point running mean of the power spectrum, result- 
ing in the plot shown in Figure 16. The top 20 peaks in the 
nighttime power spectrum are listed in Table 2. We see from 
this figure and from this table that there are two notable peaks 
in the rotational search band 10-15 year -1 , one at 12.39 year -1 , 
with power S =49.63, and the other at 1 1.27 year -1 , with power 
5=29.39. (In power-spectrum analysis, the probability of the 
null hypothesis, that the peak is due to normally distributed ran- 
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Figure 8: Gamma measurements as a function of date and time of day. The 
color-bar gives the power, 5 , of the observed signal. 



Figure 9: Temperature as a function of date and time of day. The color-bar 
gives the power, S , of the observed signal. 



dom noise, is given by e -s , where S is the power value of the 
frequency peak |9[].) 

It is interesting, and possibly significant, that the frequencies 
of the 8th and 7th peaks (11.27 year -1 and 5.54 year -1 ) are in 
close to a 2: 1 ratio. 

In previous analyses of decay data Q H 0L we have eval- 
uated the significance of a peak in the rotational search band 
by using the shuffle test 11 1 CHI . We have therefore carried out 
1,000 shuffle simulations of the nighttime data, with the result 
shown in histogram form in Figure 17. The largest value of the 
power generated in this way is 9.35, far below the actual value 
of 49.63. Hence there is no possibility that the oscillations are 
the result of independent random fluctuations in the measure- 
ments. 

It is clear from Figures 15 and 16 that the nighttime data con- 
tains oscillations that are not present, or only weakly present, in 
the daytime data, and that the difference between the daytime 
and nighttime power spectra is statistically significant. 

4. Comparison with BNL and PTB Data 

It is interesting to compare the results of our analysis 
of the GSI data with the results of similar analyses of the 
BNL (Brookhaven National Laboratory) data |Qi]] and PTB 
(Physikalisch-Technische Bundenanstalt) 11211 data. We have 
previously shown the results of simple power-spectrum analy- 
ses of these two datasets lH, H Q] . However, to facilitate visual 
comparison with our present analysis of GSI data, we now show 
instead the results of time-frequency analyses. 

We show in Figures 18 and 19 time-frequency displays of the 
BNL measurements of the decay rates of 36 C1 and 32 Si, respec- 
tively, for the frequency band 10-15 year -1 . We see that both 
of these figures exhibit features in the two bands 11.0 to 11.4 
year -1 and 12.6 to 12.9 year -1 . 

We show in Figure 20 a time-frequency display of the PTB 
226 Ra data. We see two features in the frequency range 11.0 



to 11.3 year , and another feature in the range 12.3 to 12.5 
year -1 . 

For comparison, we see in Figure 13 (formed form GSI data) 
a feature in the band 11.1 to 11.4 year -1 , and a feature in the 
band 12.2 to 12.6 year" 1 . We see that all four datasets lead to 
oscillations in the two bands 11.0 to 11.4 year -1 and 12.1 to 
12.9 year -1 . 

It should be noted that the BNL and PTB oscillations are 
intermittent in amplitude, and drift somewhat in frequency, but 
we nevertheless see broad agreement among oscillation patterns 
in all four displays. 

Bellotti et al. have recently reported the result of measure- 
ments of the activity of a I37 Cs source, as determined by an 
experiment installed deep underground in the Gran Sasso lab- 

J — L 

oratory 08J, finding no evidence of variability in that dataset. 
This result is not inconsistent with our studies to date, which 
have shown that (a) not all nuclides exhibit variability; and (b) 
among nuclides that do exhibit variability, the patterns of vari- 
ability are not all the same. 



5. Discussion 

The patterns we see in the GSI radon measurements are so 
striking as to be quite incompatible with the textbook picture 
that nuclear decay rates are all constant. The first question to be 
addressed is the extent to which the variations have an environ- 
mental or experimental origin. Figures 6, 7, 8 and 9 show that 
the gross variations of the gamma-detector measurements are 
quite similar to those of the ambient temperature. However, the 
annual waveform of the gamma-detector measurements leads 
rather than lags that of the ambient temperature. This makes 
it unlikely that the annual variation of the decay rates is due 
completely to the variation in the ambient temperature. Figure 
7 shows that the annual waveform of the gamma-detector mea- 
surements also leads that of the voltage. Hence it is also un- 
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likely that the annual variation of the decay rates is completely 
due to the voltage variation. 

On the other hand, comparison of Figures 9 and 10 shows 
that the annual variation of the gamma-detector measurements 
lags that of the solar elevation, which is compatible with the 
possibility that the annual variation of decay rates may be due 
to a direct solar influence rather than a secondary environmen- 
tal influence. Note, however, that Figure 10 represents not only 
the maximum elevation as seen at the site of the experiment, but 
also the maximum elevation as it would be seen at the diamet- 
rically opposite position on Earth. 

The frequency-hour-of-day displays in Figures 11, 12, and 
13 show that the annual oscillation is primarily a daytime ef- 
fect, but all other oscillations are primarily nighttime effects. It 
may turn out that this curious property of the time series is due 
to some experimental or environmental effect, but we have so 
far failed to identify such an effect. If this pattern proves not 
to be attributable to an experimental or environmental effect, 
indicating that it is an intrinsic property of the decay process, 
we would find it difficult to avoid the conclusion that the rele- 
vant decay process is anisotropic. There would then be a clear 
need for a dedicated experiment that could distinguish between 
isotropy and anisotropy. 

We have noted in our earlier analyses |0 01 of BNL and 
PTB data ll UN 211 that the above oscillation frequencies are too 
low to be attributable to rotation in either the convection zone 
or the radiative zone, leading to the suggestion that the solar 
core may rotate more slowly than the outer layers. This conjec- 
ture leads in turn to the suggestion that there may be an inner 
tachocline (a region involving a sharp gradient in the rotation 
rate) that separates the core from the radiative zone, analogous 
to the well known outer tachocline that separates the radiative 
zone from the convection zone. An equatorial section of the 
outer tachocline extends over the frequency range 13.7 to 14.7 
year 1 fll . 

The above similarity between the GSI oscillation frequencies 



and those found in BNL and PTB data further supports our pro- 
posal that nuclear decay rates are influenced by solar radiation 
that is somehow influenced by processes in the deep solar inte- 
rior. This leads us to look for a form of radiation that has two 
(apparently contradictory) properties: (a) the radiation can tra- 
verse the solar radiative and convection zones (and presumably 
the Earth) without being absorbed, and (b) the radiation can be 
modulated by processes in the deep solar interior and retain that 
modulation in the course of its passage through the solar inte- 
rior. 

These two requirements can be met by neutrinos: the cross 
section for scattering of neutrinos by normal matter is esti- 
mated to be extremely small, so that one would expect there 
to be negligible absorption of the neutrino flux in its propaga- 
tion from the core to Earth. However, there are believed to be 
at least three types of neutrinos involving the three flavors elec- 
tron, muon, and tau. An inhomogeneous magnetic field (and 
magnetic fields are necessarily inhomogeneous) can lead to fla- 
vor oscillations by the RSFP (Resonant Spin Flavor Precession) 
process ll4lll5Ul6ll . One can imagine that such a modulation 
might be imprinted on the neutrino flux in the deep solar inte- 
rior and remain imprinted in the flux during its passage through 
the outer layers of the Sun. 

The obvious and crucial objection to this proposal is that 
there is no known mechanism by which neutrinos can influence 
beta decays. Further experiments and future theoretical work 
will hopefully lead to a resolution of this paradox. 
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Hour of Day 



Figure 12: Frequency-Hour-Of-Day display for the frequency range 5-10 
year . The color-bar gives the power, 5 , of the observed signal. 

Appendix A. Experimental Setup 0] 

The experimental setup (Figure 21), located in the yard of 
the Geological Survey of Israel (GSI) Laboratory in Jerusalem, 
is contained in a square tank comprised of welded 3-mm-thick 
iron plates. The tank is sealed from the outer environment. The 
upper part of the tank contains air and the lower part contains 
ground phosphorite. Uranium ( 238 U) in the phosphorite pro- 
duces radium ( 226 Ra) that in turn produces radon ( 222 Rn). A 
gamma-ray sensor (gamma-C), located in the upper part of the 
tank, is contained in a lead pipe (length 30 cm, outer diame- 
ter 10.5 cm., thickness 0.5 cm), which has a perforated lead 
plate at its lower end. This plate serves to reduce the direct 
gamma radiation from the phosphorite, while allowing gas to 
enter the gamma-detector chamber. The pressure and temper- 
ature sensors, a datalogger, and a 12-volt battery are located 
in a box under the shack. The gamma detector is a 2x2 Nal 
detector (PM-11 detector, ROTEM Inc., Israel), which records 
gamma-ray impulses of energy above 50 keV. Measurements 
are recorded at 15-minute intervals. 

Appendix B. Environmental Drivers of Radon Production 

The variation patterns of radon in geogas (subsurface air) 
have been observed and investigated since the second half of 
the last century. Numerous publications aim to clarify possi- 
ble reasons for observed radiation patterns and their potential 
geodynamical significance. Due to its noble-gas character, its 
constant generation, and its radioactive decay, radon provides a 
unique indicator of subtle variable processes in the natural en- 
vironment. It is significant for geological processes that radon 
often shows large and complex variations that have no obvi- 
ous atmospheric causes. These variations are believed to reflect 
subsurface processes. However, there is as yet no comprehen- 
sive theoretical understanding of observed radiation patterns. 
This lack of understanding of the causes of radon variability 
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Hour of Day 

Figure 13: Frequency-Hour-Of-Day display for the frequency range 10-15 
year" 1 . The color-bar gives the power, 5 , of the observed signal. 

hampers its utilization as a proxy for active geodynamics. For 
this reason, it is important to try to understand radon variability 
that is not obviously related to atmospheric variability. 
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Figure 14: Power spectra over the frequency range 0-5 year -1 . Power from Figure 15: Power spectra over the frequency range 5-20 year - '. Power from 
daytime data shown negative, and power from nighttime data shown positive. daytime data shown in blue, and power from nighttime data shown in red. 



Table 2: Top 20 peaks in the 41 -point running mean of the power spectrum 



formed from nighttime data. 


Frequency (year 1 ) 


Power 


Order 


1.99 


286.22 


1 


0.37 


137.00 


2 


1.06 


130.99 


3 


2.89 


57.21 


4 


12.39 


49.63 


5 


5.00 


35.65 


6 


5.54 


30.27 


7 


11.27 


29.39 


8 


8.90 


29.38 


9 


17.54 


22.18 


10 


10.76 


20.12 


11 


7.96 


17.87 


12 


7.80 


17.28 


13 


4.32 


16.46 


14 


18.53 


16.07 


15 


12.98 


15.52 


16 


13.62 


15.34 


17 


9.80 


15.01 


18 


18.15 


15.00 


19 


19.68 


13.41 


20 



60 



40- 




Frequency (year 1 ) 



Figure 16: 41-point running mean of the power spectra over the frequency range 
5-20 year -1 . Power from daytime data shown in blue, and power from night- 
time data shown in red. 
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Figure 17: Histogram formed from 1,000 shuffle simulations of the 41-point 
running mean of the power spectrum formed from nighttime data. None has 
power in the rotational band as large as 49.63. The largest value in the simulated 
data is 9.35. 
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Figure 18: Frequency-time display of the CI decay measurements made by 
the BNL (Brookhaven National Laboratory) experiment. The color-bar gives 
the power, S , of the observed signal. 
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Figure 19: Frequency-time display of the 32 Si decay measurements made by 
the BNL (Brookhaven National Laboratory) experiment. The color-bar gives 
the power, S , of the observed signal. 
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Figure 20: Frequency-time display of the 226 Ra decay measurements made by 
the PTB (Physikalisch-Technische Bundesanstalt) experiment. The color-bar 
gives the power, S , of the observed signal. 
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Figure A. 21: Experimental layout, showing the locations in the tank of the 
phosphorite (the source of 222 Rn) and the radiation detectors. This article ana- 
lyzes measurements made with the gamma-C detector, which is housed located 
in a lead tube. Pressure and temperature sensors, a data logger, and a continu- 
ously charged 12 V battery are located in a separate enclosure attached to the 
wall of the shack. 
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